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Abstract. We investigate conditions for logarithmic complete mono¬ 
tonicity of a quotient of two products of gamma functions, where the 
argument of each gamma function has different scaling factor. We give 
necessary and sufficient conditions in terms of nonnegativity of some 
elementary function and more practical sufficient conditions in terms 
of parameters. Further, we study the representing measure in Bern¬ 
stein’s theorem for both equal and non-equal scaling factors. This leads 
to conditions on parameters under which Meijer’s G- function or Fox’s 
//-function represents an infinitely divisible probability distribution on 
positive half-line. Moreover, we present new integral equations for both 
G-function and //-function. The results of the paper generalize those 
due to Ismail (with Bustoz, Muldoon and Grinshpan) and Alzer who 
considered previously the case of unit scaling factors. 
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1. Introduction. Recall that a nonnegative function / defined on (0, oo) is called com¬ 
pletely monotonic (c.m.) if it has derivatives of all orders and (— l) n /i n \x ) > 0 for n > 1 
and x > 0 [20, Definition 1.3]. This inequality is known to be strict unless / is a constant. 
By the celebrated Bernstein theorem a function is completely monotonic if and only if it is 
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the Laplace transform of a nonnegative measure [20, Theorem 1.4]. The above definition 
implies the following equivalences 


/ is c.m. on ( 0 , oo) yy / > 0 and —f' is c.m. on ( 0 , oo) 

yy —f is c.m. on ( 0 , oo) and lim f(x) > 0 . ( 1 ) 

x —>00 

A positive function / is said to be logarithmically completely monotonic (l.c.m.) if 
— (log/)' is completely monotonic [20, Definition 5.8]. According to (1) 


/ is l.c.m. on ( 0 , oo) yy (—log f(x)) r > 0 and (log f)" is c.m. on ( 0 , oo) 

yy (log/)" is c.m. on ( 0 , oo) and lim (—log > 0 . ( 2 ) 

£—>•00 

The class of l.c.m. functions is a proper subset of the class of c.m. functions. Their 
importance stems from the fact that they represent Laplace transforms of infinitely divisible 
probability distributions, see [20, Theorem 5.9] and [19, Section 51]. 

The study of complete monotonicity of the ratio 


u(x)= n 

i= 1 


T(a; + a*) 

r {x + biY 


( 3 ) 


where T stands for Euler’s gamma function and p > 1, has been initiated by Bustoz and 
Ismail who demonstrated in their 1986 paper [4] that for p = 2, a± = 0 and <q + a 2 = b\ + 62 
this function is l.c.m. on (0, 00 ). Eight years later Ismail and Muldoon showed in [7] that 
U(x ) is l.c.m. on (0, 00 ) for general p if Y^=\ a i = Sf=i an( l — b 2 = ■ ■ ■ = b p > 0, 
> 0, or b\ = b 2 = ■ ■ ■ = b p _ 1 = 0, b p > 0, a* > 0. Their result is, in fact, formulated for 
the ratio of (/-gamma functions, but the proof works for U(x) just as well. The subject was 
further pursued by Alzer [ 1 ] who showed in 1997 that U(x) is l.c.m. on (0, 00 ) if 


0 < ai < a ,2 < • • • < a p , 0 < bi < b 2 < ■ ■ ■ < b p , 

k k 

and a % < bi for k — 1, 2 ..., p. 

i =1 i— 1 


( 4 ) 


These inequalities are known as weak supermajorization [14, Definition A.2] and are abbre¬ 
viated as b -< 11 a, where a= (ai,..., a p ), b = (b\,... , b p ). In their 2006 paper [5] Grinshpan 
and Ismail found new sufficient conditions for U(x) to be l.c.m. of when p = 2 n or p = n\/ 2, 
n > 1. We will explain and slightly generalize their results in the next section. Finally, 
in 2009 Guo and Qi [6] used another approach to investigate logarithmic complete mono¬ 
tonicity of U(x) for arbitrary real values of a ,, bi. Their results, however, lead to complete 
monotonicity of U(x) on some subinterval of (0, 00 ) of the form (7, 00 ), where 7 > 0. 

The purpose of the present paper is to investigate complete monotonicity of the ratio 


W{x) 


nr=i r (a^ + ai) 
nj =1 r (BjX + bjY 


( 5 ) 


where A = (Ai,...,A p ) and B = (B \,..., B q ) are strictly positive scaling factors, while 
a = (d!,..., a p ) and b = ( 6 1; ..., b q ) are nonnegative. We find a necessary and sufficient 
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condition for W to be l.c.m. on (0, oo) in terms of nonnegativity of some function which, 
unfortunately, is not easy to verify. Further, we supply several simple sufficient conditions 
for such nonnegativity in terms of the vectors A, B, a, b as well as some necessary conditions. 
When W is completely monotonic, we proceed by deriving its representing measure in the 
Bernstein theorem which leads us to studying some new properties of Fox’s //-function. We 
begin, however, by revisiting the ratio defined in (3) which we call the unweighted case. In 
the following section we refine some of the known results for U(x) and discuss its representing 
measure in Bernstein’s theorem. 


2. The unweighted case revisited. All proofs of complete monotonicity of U(x) we are 
aware of are based on the integral representation of the logarithmic derivative of the gamma 
function [13, p.16] 


'ijj(z) = 


n 


F(*) 


= -7 + 


-°o e -t 


— tz 


1 - e 


— t 


-dt, 


( 6 ) 


where 7 is Euler-Mascheroni constant. The next two properties of digamma function -0 will 
be repeatedly used in the sequel 


$( x ) = I YZ e -t dt ’ x > °> ( 7 ) 

0 

and 

l 

i>(z) = log(z) - — + 0(z~ 2 ) as z ->■ 00, (8) 

see [2, Theorem 1.6.1, Corollary 1.4.5]. 

Representation (6) lead Grinshpan and Ismail to formulate Lemma 2.1 in [5] stating that 
q(x) = nLi + Afc) with Y^k=i = 0 and A*, > 0 is l.c.m. on (0, 00) if and only if 

p 

v(t) = ^2 a kt Xk >0 for t G (0,1]. 
fc=i 

In fact, looking at asymptotic expansion of —(logq , (a;)) , as x — > 00 it is easy to sharpen this 
lemma as follows: q[x) is l.c.m. on (0, 00) iff v(t) > 0 on (0,1] and Ylk=i a k = 0. When 
q(x) = U(x), where U is defined in (3), we have 

p 

U{x) is l.c.m. on (0, 00) v(t) = ^2(t ak ~t bk )> 0 on (0,1], (9) 

k =1 

As mentioned above Alzer noticed that v(t) > 0 on (0,1] under majorization conditions 
(4) which follows from 1949 result of Tomic [14, Proposition 4.B.2], Let us remark that 
for p = 2 conditions 0 < min(ai,a2) < min(6i,6 2 ) and a\ + a 2 < 5i + b- 2 , equivalent to 
(4), are necessary and sufficient for nonnegativity of v(t). Necessity follows by considering 
the asymptotics of v(t) as t —>■ 0 and t —> 1. A proof of a more general result is given in 
Corollary 1 below. 
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Grinshpan and Ismail took another path in [5] and considered two types of factorization 

n 

v(t) = (t aj - t ai ) and v(t) = JJ(1 - t ai ). 

i —1 

It clear that for a\ > «2 > • • • > ot n > 0 for the hrst factorization and a* > 0 for the 
second we get expressions nonnegative on (0,1]. The corresponding value of p is n\/ 2 for the 
hrst factorization and 2 n ~ 1 for the second. In both cases the authors of [5] found explicit 
combinatorial descriptions of the vectors a, b in (3) which lead to the above factorizations. 
We would like to remark here that taking v in the form 


v(t) = ]J(t 0i - t ai ) 

i= 1 


with a* > fa > 0 for i = 1,..., n certainly results in nonnegative v on (0,1]. This, of course, 
can be reduced to the second factorization of Grinshpan and Ismail by factoring out 
Also one can get the hrst type of factorization by properly choosing fa. The vector a (b) is 
expressed in [5] in terms of sum over even (odd) permutations of the numbers a, for the hrst 
factorization, and in terms of sums of the numbers a* indexed by components of increasing 
integer vectors for the second factorization. In the theorem below which generalizes both 
factorizations, the description of the vectors a and b in terms of the numbers and fa turns 
out to be different and simpler. 


Theorem 1 Let I = {1,2, and suppose X odd ( X even ) comprises all subsets of I 

having odd (even) number of elements, 0 G T even . Suppose that a t > fa > 0 for i — 1 ,..., n. 
Then 

Yljeleven ^ ( X + + J2i£l n \j A 

U(x) = 


n 


J&lodd ■ 


r ®+E <e ja* + E 


'ielnV A 


is l.c.m. on (0, oo). 


Proof. We will write / = I n , T odd = I odd (n), T even = I even (n) to emphasize the depen¬ 
dence on n. According to (9) it suffices to verify that v(t) > 0 on (0,1], where 

. ^ E E hi . ^ E «;+ E hi 

Jel even (n) J&l 0 dd(n) 


We will demonstrate that 

n 

v(t) = UT‘ - n 

i= 1 

by induction in n. Indeed, the case n — 1 is obvious. The induction step now follows from 
the identity 



E a i + E hi 

i£ln\J 


J EX e , 


n(n) 



E X Pi+Pn 



X Cx.i+Ot.n-\- XI Pi 


Jel even (n-l) 


^eX odd (n— 1) 
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and similar identity for the negative term in v(t). Nonnegativity is now obvious from cq > 
Pi > 0 for i — 1 ,..., n. □ 

The next example generalizes [5, Corollary 1.4] 

Example 1. For cq > Pi > 0 for i = 1, 2, 3 the function 


r(x + Pi + P2 + P^T^x + Pi + ot.2 + a 3 )r(a; + ai + P2 + o: 3 )r(:r + oq + 02 T / 3 3 ) 
r(x + a.\ + P2 + / 3 3 )T(x + Pi + «2 + / 0 3 )r(x + /?i + P2 + o: 3 )r(:r + oq + o;2 + cx 3 ) 


is l.c.m. on (0, oo). 

Further, changing variable in [18, formula 2.24.2.1] or in particular case of [11, Theo¬ 
rem 2.2] we get 


n r(x + ap 
_ T(x + bi) 



(0,oo) 



( 10 ) 


provided that YP!= i (h ~ ap 
integral 


^p,p 


> 0. Here G p p f p is Meijer’s G'-function defined by the contour 

b\ = J_ P T(ai + s)... r(a p + s) 

ay 2ni J T(bi+s)... T(b p +s) 
c 


This function is a particular case of Fox’s //-function defined in (21) below on setting p = q 
and Ai = Bi = 1, % = 1,... ,p in that definition. See text below (21) for a description of the 
contour C. It follows from an expansion due to Nprlurid [16, (2.28)] that for YPPi=\ P>i~ a P — 0 
formula (10) must be modified as 


TT r(x + ap 

fjr (x + bp 


= / {<^0 + ( e “* 

[0,oo) 



(ii) 


where 5 o denotes the unit mass concentrated at zero. More information about Nprlund’s 
expansions can be found in our forthcoming paper [9]. Representation (10) was previously 
observed by us in [10]. 

The l.c.m =>■ c.m. implication mentioned above U translates into the next assertion: 


v(t) =^2(t ak - t bk ) > 0 on [0,1) 


/-ip, o 


k=1 


X 


> 0 on (0,1) 


under additional conditions a* > 0, i — 1 ,... ,p, and ~ a P > 0- The following much 

stronger assertion is supported by numerical evidence: 


Conjecture 1 If a* > 0, i — 1,... ,p and YPi=i (h — ap > 0 then 


#<i zeros of G p p ’° p ( x 


on (0,1) > < #< zeros of v(t) on (0,1) 


We think that this conjecture is encoded in the following integral equation for Meijer’s G- 
function which we believe to be new. 
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Theorem 2 Suppose a, b > 0 and fi > 0. Then Meijer’s G-function satisfies the follow¬ 
ing integral equation: 


\og{l/x)G 


p, o 

p,p 



x bk \ dt 
t bk ) t — x 


for 0 < x < 1 . 

We postpone the proof of this result until proof of Theorem 8 , containing a more general 
statement. Using (2) we were also able to demonstrate that if G^f p (x) has a zero on (0,1) 
then v(t) also has at least one zero on this interval. 


3. Generalized gamma ratio. For the positive vectors A = (A l5 ..., A p ). B = (U 1; ..., B q ) 
and nonnegative vectors a = (ai,..., a p ), b = (bi,... ,b q ), consider the positive function 


n:.i r i-i,. + at 

= Ii r (B jX + b s ) 


defined on (0, oo). Our goal is to study logarithmic complete monotonicity of W. It is 
convenient to start with 


Lemma 1 The function (log W)" is c.m. if and only if 

P p—diu/Ai Q p—biii/Bi 

P ( U ) = Y 1 _ e -u/ Ai - Y I _ e -n/B t ^ 0 f° r dl u > °- ( 12 ) 

i =1 i =1 

In the affirmative case 

OO 

(log W)" = J e~ xu uP(u)du. (13) 

o 

Proof. Differentiating (6) we have 

p Q 

(log Wf = Y, AWiAx + afi~Y + bfi. 

i =1 j =1 

Using (7) and making change of variable u = Afi or u = Bfi in the appropriate integrals we 
get 


(log W)" = Y 1 A] 


"°° te~( AiX+ai ' )t dt 


i =1 


1 — e 


-t 


j =i 


f°o te -{ Bi x+b s yt dt 




ue -xu-u ai /A idu 9 

_ p u/Ai / ^ 


1 — e t 

50 ue- xu ~ ub ^ B ^du 


e~ xu u 


P ^ e ~aiu/Ai 

t ^ 1 _ P u/A-i 


1 - e~ u / B . 
e -bju/Bj 


3 = 1 J0 


i— 1 


5Z ! z e -u/Bj \ du = i e XU uP(u)du. 


3=1 
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The integral converges around zero due to asymptotic formula (19) below. Convergence 
at infinity is obvious. Hence, (log W)" is c.m. iff P(u) > 0 by Bernstein’s theorem [20, 
Theorem 1.4]. □ 

Remark 1. By substituting t = e~ u condition (12) can be also written in the form 


Q(t) = J2 


2—1 


fCLi/Ai 

1 — t x ! Ai 



fbj/Bj 

1 — 


> 0 


for t G ( 0 , 1 ). 


(14) 


If (log W)" is completely monotonic two interesting possibilities appear: (log IT)' > 0 
and (log IT)' < 0. These two possibilities are handled in the next two theorems. Recall that 
a nonnegative function / on (0, oo) is called Bernstein function if /' is c.m. [20, Definition 3.1 
and below]. 


Theorem 3 Suppose A, B > 0 and a, b > 0. Then the function 

p <2 

(log IT)' = 22 Af{Ax + a,') - 22 + bj) 

i=l j =1 

is a Bernstein function iff condition (12) holds and 

i^o {El =1 A ^) - Y ? j= 1 > o. 

Proof. The proof is immediate from Lemma 1 and the following equivalences 


/ is Bernstein function f is c.m. and / > 0 on (0, oo) 

f is c.m. on (0, oo) and lim fix) >0. □ (15) 

x^O 


Theorem 4 The function W is l.c.m. if and only if 


= 11 II' ; 


and condition (12) holds. In the affirmative case 


-(log W)' = I 
Jo 


e ~ xu p[u)du + log(l/p). 


(16) 


(V) 


Proof. We have 

q p 

-(log IT)' = ^2 B j if(B j x + bj) - 22 AfiAx + af). 

j = 1 7=1 


Further it follows from ( 8 ) that 


f(Cz + c) 


log 2 + log C + 



1 

Cz 


+ 0{z~ 2 ) 


as z —)■ oo. 
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Substituting the latter formula into the former, we get 


-(log WO' 


iog(z) \ A > 

. j =1 i =1 


E B > m^) - E A log( A 

.3 = 1 


i= 1 


l 


- Y1 ai + 

. t=i ' i 


p-g 


+ 0(x 2 ) as x —» oo. (18) 


If conditions (12) and (16) hold than 

lim (— log W(x))' = — log p > 0. 
£—>•00 


Then Lemma 1 and relations (2) imply that W is l.c.m. In the opposite direction, if W 
is l.c.m. than — (log IT)' is c.m. and hence decreasing which in view of (18) implies that 
EL, A > On W ie other hand, relations (2) show that lim 3 ,_ >00 (— logW(x))' > 0 

so that Ef =i Ai < Ei=i A again by (18)- This proves the first condition in (16). Once this 
condition is satisfied nonnegativity of the limit leads to the second condition in (16). Finally, 
integrating (13) we get 


(loghF(x))'= f (log W(x))"dx + log p 

J OO 


oo 

J uP(u)du 

o 


e tu dt + log p 


OO 


OO 

o 


log(l/p). 


□ 


Remark 2. Expression of the form S\ = 'YJ ‘ = , A, log(A) for positive numbers A is 
known as Shannon’s entropy in information theory, so that the second condition in (16) can 
be restated as ’’entropy of the vector A does not exceed that of the vector B”. 

We now collect the necessary conditions for logarithmic complete monotonicity of W in 
the next corollary. 


Corollary 1 The following conditions are necessary for W to he logarithmically com¬ 
pletely monotone: 


a) = 

b) p < 1 

c) t = E?=i b J ~ E?= i a i + \(p- q)> o 

d) min (m/A) < min (bi/BA 

1 <i<p 



Proof. Necessity of a) and b) has been demonstrated in Theorem 4. Next, straightfor¬ 
ward computation shows that 

1 f p i i p i ) 

p( u ) = - < a - B i f b i ~ ai + 2 ^ ~ i + as u ° - ( 19 ) 

u l i= 1 3=1 J {3=1 i=1 J 

Since the first sum is zero by a) we conclude that c) is necessary for nonnegativity of P(u). 
Further, using representation (14), we easily compute 

Q(t) = t a (l + o(l)) - t p (l + o(l)) as t —> 0, 


where a = min ( at/Ai ) and (3 = min ( bj/Bj ). Hence, we conclude that condition d) is also 

l<i<p 1 <j<q 

necessary. □ 

Remark 3. Another method to show the necessity of condition p, > 0 (condition c) 
above) is as follows. Calculating further derivatives of / and using 

= AzA + _A_ + 0(x— 2 ) as x —f oo, 

we get 


{-l) n f {n \ X ) = (-l) n+1 ^(4^V (n) (As + ai) - (-l) n+1 B] +1 ^ n \B jX + bj)) 

3 = 1 


i= 1 

f- f («- mri i n! - 4 - +1 

(AiX + ai) n 2 (AiX+di) 


n+1 


-E 


(n - 1)!R” +1 


, + n'-Br 1 


1 (BjX + bj) n 2(B jX + bj) 


n+l 


>+0(x 


—n—2\ 


E 

i =1 


(n - 1)!A 


x 7i 


na, \ n\ 
1-^1 + 


(n — 1)! 


AiX J 2x n+1 
q 


-E 

j=i 


(n — l)!Rj 


x 71 


1 - 


nb-i 


+ 


n! 


BjX / 2 x n+1 


+0(a; 


—n—2\ 


X n 


^ _ B i f + 
*=i i=i 


n! 


x 


n+l 


b i ~ S + o “ ?) f + °( a; " 


—n—2\ 


TV. 


,j =1 *=1 

9 P 




■n+l 


X) b f ai + ~ q ) }+°( 

. 3 =1 i=1 




—n—2\ 


where the first sum vanishes due to the necessary condition a) from Corollary 1. This shows 
that this must require that Ej=i ~ EL i a * + ^(p ~ Q) — 0 m order for (—1 ) n f( n \x) to be 
nonnegative for large x. 

Given the vectors A, B, a, b condition (12) is not easy to verify. The next theorem 
provides some practical sufficient conditions. 

Theorem 5 Inequality (12) is true if any of the following conditions holds: 


(a) E )=i B j = ELi A % and max {a l /A i ) < min (b 3 - 1 )/B j - 

1 < 2 <P 

(b) p = q, ELi B i = Ef= i mth Ai> Bi for i = 1,... ,p - 1, and 

max b k /B k < (b p - 1 )/B p , a^Ai < (h, - 1 )fB i for i = 1,... ,p; 

l<k<p—l 
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(c) p = q, 0 < a-i/Ai < CL 2 /A 2 < ... < a p /A p , 0 < bi/Bi < b 2 /B 2 < ...< b P /B p , 


0 < 1 / A\ < 1/A 2 < ... < 1/Ap, 0 < 1/Bi < 1/B 2 <■■■< 1/B P and 


Yh= 1 a i/A < Y!l=] bi/Bi, Et, 1/A < Ell 1/A /or eac/i integer k = 1 ,... ,p. 

Proof, (a) Denote a' = a t /A v , A[ = 1/A, b[ = bi/Bi , A = 1 /Pj, x = e u . Then by the 
mean value theorem we have 


P e-a iu/M P 

/ ^ _ g w/^ 


1 - 


v_-_= V 

x a\ _ 


E 4 - 


2=1 2=1 2—1 
where q G (a' — A, a D- Hence, for P(tt) we get 

1 1 P 


—' A[x Ci log x log x 




p (A = - 


—dnU 


( 20 ) 


2=1 


2=1 


where dj is some point of the interval (&' — B'-, P'). Denote c = maxq, <2 = min d r Suppose 
max(aj/A) < min( 6 j — 1 )/Bj, which is equivalent to max(a') < min( 6 ' — B'-). Then c < d 
and nonnegativity of P(u) follows from the relations 


Y A iX ~ Ci > Y A i x ~° > Y A i x ~ d = Y B i x ~ d 


x 


-dj 


i= 1 


i =1 


i= 1 


2=1 


2=1 


p p 

(b) Denote S=E5j = El' Nonnegativity of P(u) follows from representation (20) 

2 =i *=1 

and the following chain 

p p p -1 / p-i \ 

Y A ie - CiU > Y Ae~ diU = Y A e~ diU + \ S ~ Y Ai ) e ~ dpU 

2=1 2=1 2=1 \ 2=1 / 

p —1 p —1 p 

= Y A(e“ diU - e~ dpU ) + P^e - *" - e _dpU ) + S , e _dpU = Y Ae“ di E 

2=1 2=1 2=1 

The first inequality here is due to conditions cq/A < (b—1 )/Pj, i = 1 ,..., p which imply that 
q < di. The second inequality follows from conditions Ai > B, , i = 1,... ,p — 1, combined 
with inequality max b^/B^ < (b p — 1)/P p which ensures that each term e~ diU — e~ dpU is 

i < fc < p—1 

nonnegative. 

(c) Expanding each expression like (1 — e _n / Ai )^ 1 in geometric series and exchanging the 
order of summations we get 


OO P 


p ( m ) = ]eS( g— aju/Ai-lu/Ai _ ^-bju/Bi-lu/Bi j 


3 — aju/Aj-lu/Ai _ —bju/Bi—lu/Bi 
V 

1=0 i= 1 

Conditions c) guarantee that for any l > 0, u > 0 we have 

( b\u/B\ + lu/Bi,b p u/B p + lu/B p ) -A (oq-u/A + lu/Ai,a p u/A p + lu/A p ) 
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(see (4) for definition of weak supermajorization -< w ). Hence, according to [14, Proposi¬ 
tion 4.B.2] we conclude that 

p 

- aiu/Ai—lu/Ai _ ^-biu/Bi-lu/Bi \ ^ q 

i= 1 



since x —> e~ x is convex and decreasing. □ 

Remark 4. Conditions (c) from Theorem 5 are only compatible with condition i A? = 
Bj if A = B (up to permutation). Indeed, majorization B' -<" A (recall that 
A'- = Aj 1 , Bj = Bj l ) forms a part of conditions (c). According to [14, Proposition 4.B.2] 
this majorization implies that 


Y. b > 


3 =1 





3 = 1 


because the function x —> 1/x is decreasing and convex. Further, [14, 3.A. 6 .a] says that 
for continuous strictly decreasing functions the above inequality is strict unless A is a per¬ 
mutation of B'. This brings us to the conclusion that for A ^ B (modulo permutations) 
conditions (c) can be used to check whether (log Wy is a Bernstein function using Theo¬ 
rem 3, but they cannot be used to check whether W is completely monotone. On the other 
hand, if A = B, conditions (c) reduce to checking the majorization a' -< w b' (a)- = ajAj 1 , 
bj = bjAj 1 ). This majorization is weaker than (a) or (b) in this particular situation. One 
reason why the case A = B might be important is applications in probability as explained 
in Remark 5 below. 


Example 2. The function 

t r( 2 a; + 0.4)r(3a; + 2A)T(x + 0.9) 

X ^ r(a: + 2)r(5a; + 6) 

is logarithmically completely monotone since A\ + A 2 + A 3 = B 1 + B 2 , p = 0.03456 and 
max 1 <j< 3 (aj/Aj) = 0.9, min 1 <j< 2 (( 6 j — 1 )/£>,:) = 1. Hence, necessary conditions (16) are 
satisfied and P(u) > 0 according to Theorem 5(a). Note, that neither conditions (b) nor (c) 
from Theorem 5 can be applied here as p 7 ^ q. 

Example 3. The function 

v r( 2 x + 0.8)r(3a: + 8)T(1.4x + 2.3) 

1 ^ r(x + 1.5)r(2.4x + 7.8)r(3a: + 11) 

is logarithmically completely monotone since Ai + A 2 + A 3 = B\+B 2 + B 3 , p = 0.783668, so 
that necessary conditions (16) are satisfied and A 2 > Bi, A 2 > B 2 , maxi<j< 2 bi/B^ = 3.25 < 
(63 - 1 )/B 3 = 3.33, a l /A 1 = 0.4 < (b x - 1 )/B 1 = 0.5, a 2 /A 2 = 2.66 < (b 2 - 1 )/B 2 = 2.83, 
a 3 /A 3 = 1.643 < ( b 3 — 1 )/B 3 = 3.33, so that P(u) > 0 by Theorem 5(b). Note that 
conditions (a) and (c) of Theorem 5 are violated here. 

Example 4. The function 

x —> 4y>(4x + 0.7) + 2ijj(2x + 1.8) — 3^(3x + 0.6) — ij}{x + 1.2) 
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is a Bernstein function, since A\ = 4 > A 2 = 2, B\ = 3 > B 2 = 1, A 1 1 = 0.25 < B l 1 = 0.33, 
Af 1 + Af 1 = 0.75 < Bf 1 + Bf 1 = 1.33, aiAf 1 = 0.175 < b x Bf 1 = 0.2, a x A 7 1 + a 2 Af 1 = 
1.075 < b\Bf l + b 2 Bf 1 = 1.4, so that P(w) > 0 by Theorem 5(c) and the claim follows 
from Corollary ??. Note, that the function W with parameters from this example is not 
logarithmically completely monotone since necessary conditions (16) are violated. 

Example 5. The function 

t T(3x + 0.8)T(2.2a: + 1.8)T(1.4x + 2.3) 

:1 ^ T(3x + 1.2)T(2.2x + 1.7)T(1.4a; + 2.5) 

is logarithmically completely monotone by Theorem 5(c). Indeed, Ai = Bj, i = 1,2,3, and 
majorization conditions are easily verified. 


4. The representing measure. 

defined by 

(B,b) 
(A, a) 


H p ’° 

Q,P 


We will need a particular case of Fox’s //-function 
nLi r (A k s + a k ) s 


1 

2iri 


n?.,r(B iS + 6 i ) 


■z S ds, 


( 21 ) 


where A k , Bj > 0 and a kl bj are real. The contour C can be either the left loop £_ starting 
at — 00 + ia and ending at — 00 + i/3 for some a < 0 < /3 such that all poles of the integrand 


lie inside the loop; or the right loop C + starting at oo + ia and ending at 00 +i/3 and leaving 
all poles on the left; or the vertical line C ic , s Rz = c, traversed upwards and leaving all poles 
of the integrand on the left. Denote the rightmost pole of the integrand by 7 : 

7 = - min ( a k /A k ). 

l</c<p 

Recall the definition of p from (16) and the definition of // from Corollary 1(c): 


/- = n c n D 

k =1 j =1 


3 


<1 


d = 2^ b i 
3 = 1 


ira k + P - 

k= 1 


Existence conditions of //-function under each choice of the contour C have been thoroughly 
considered in the book [11], Under restrictions (16) and x > 0 Theorem 1.1 from [11] states 
that H£p(x) exists if we choose C = C + or C = £j C under additional restriction p > 1. Only 
the second choice of the contour ensures the existence of the Mellin transform of H^{x) 
as demonstrated in [11, Theorem 2.2], In our next theorem we relax the condition p > 1 
to p > 0 and demonstrate that the first condition in (16) leads to the compact support of 

HISix). 

Theorem 6 Suppose p > 0 and Y^=\ A = 1 Bj. Then the integral over the contour 

Ci c with c> 7 in the definition of H//^(x) converges for x > 0 except possibly for x = p and 


H'Pfi 

q,p 


x 


(B,b) 
(A, a) 


= 0 


( 22 ) 


for x > p. Moreover, under these restrictions, the Mellin transform exits for > 7 and 


H p '° 

q,p 


x 


(B,b) 
(A, a) 


x 


s— 1 


dx = 


nLi r(AfcS + Ofc) 
npWTy' 
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Proof. By applying Stirling’s asymptotic formula and in view of A k — Sj=i Bj it is 

not difficult to derive the formula [17, (2.2.4)] 


m=i r(^4fc<s + a*;) 

HL%Ty 


= Ap s s^ + p s g(s) : 


with g(s) — 0(s M x ) as |s| —* oo uniformly in | arg s\ < tt — e for any 0 < e < 7r and 

A = (2n) ip - q)/2 Y\ P Al k ~ 1/2 T\ q Bt /2 ~ bj . 

This asymptotic behavior implies that t —> g(c + it) is absolutely integrable continuous 
function on the real line for any c > 7 so that the integral v(x) = Ic~i^( x /P) S 9(s)ds 
exists and we are in the position to apply the Mellin inversion theorem yielding 


Denote 


x s 1 v(x)dx = p s g(s). 


h(x) = — / Ap s s~ fl x~ s ds. 


An application of [3, ( 6 ), §12] after a simple rearrangement gives 


h(x) = 


+OO 

A f glog (p/ x )(c+it) 

27 T J (c + itY 


dt — < r(/i) 


(!ogf) M 1 0 < x < p, 


0 , x > p. 


It then follows from (23) that 


x 


(B, b) 
(A, a) 


_L / nLl + gfe) s d Jl(x') + v(x) 
2m J UU I '(B j s + b J ) X 


and the above integral exists for all x > 0, x 7 ^ p. Direct computation gives 


x^H^ 


(B h) W — 

(A, a) dX Tip 


PY~ l 1 


x s dx + p s g(s ) = Ap s s M + p s g(s), 


which is equivalent to 


x‘-'hz:°(x 


(b, b) \ dx = n;;.i r(/i tS +o. t > 


(A, a) 


n?=. r(B jS + bj 


Next we show that (22) holds so that the integration in the above formula is in fact over the 
interval (0, p). To this end we apply an argument similar to Jordan’s lemma. According to 
the definition ( 21 ) 


(B, b) 


= lim 


c-\-iR 

i f rifc=i r(^AA + Ofc) _ 


(A, a) J r^-oo 2iri J n?=i r(B?s + bj 


z S ds. 
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By Cauchy’s theorem the last integral equals to the integral along the right semicircle of 
radius R centered at c. Hence, we need to prove that 


7r/2 



-tt/2 


i r(j4feC + AkRe 1 ^ + <ifc) 

n ’ =1 r (BjC + B s Re+ bj 


g — (c+R cos (p+iR sin if) log x +iip^^ 


goes to zero as R —* oo for x > p. Setting s = Re llf in (23) we get 

mV 1 T ^ A a kC + + = Ap^iRe^y^l + OiR- 1 )) as R ^ oo. 

rij =1 r(.BjC + BjRe l<fi + bj) y J y y ” 

Hence, for sufficiently large R and some constant M\ 

t/2 


1^.1 < 


M x e 


— clog# 


2 vr 


R 


-/X+1 cosip 




—7r/2 


Employing the elementary inequality cos 99 > 1 — -<p, 0 < <p < tt/2, we obtain the following 
estimate with some constant M 2 (recall that x > p): 


t/2 


\I R \ < 2M 2 R- fi+1 / e - R{1 ~^ )log - P dip = nM 2 R~ fl 


- u ( 1 ~ e 


-R log; 


logf 


Hence, lim Ir — 0 which completes the proof of the theorem. 

R—> oo 

Theorem 7 Suppose that p > 0 and 

w( X ) = 


is l.c.m. Then 


W(x) = 


n; =1 r {Bjx+bi) 

B, b) 


e~ tx HP§ ( e - 4 


(A, a) 


dt, 


□ 


«/log(l /p) 

and the H-function in the integrand is nonnegative. In particular, the conclusion is true if 
conditions (16) and hypotheses of Theorem 5 are satisfied. 


Proof. According to Theorem 4 logarithmic complete monotonicity of W(x) implies 
Y^k= 1 A-k = Bj, so that we are in the position to apply Theorem 6 yielding the formula 


W{x) 


ELL r (AiX + af) 

n; =1 r (Sjx+bj) 



(B, b) 
(A, a) 


u x 1 du. 


The claimed Laplace transform representation for W(x) follows by substitution u = e“L 
Nonnegativity of HJf^(u) on (0,p) follows from Bernstein theorem in view of uniqueness of 
the measure with given Laplace transform, see [20, Proposition 1.2] or [22, Theorem 6.3]. □ 
Theorem 7 requires /i to be strictly positive, while the necessary conditions from Corol¬ 
lary 1 allow for pL = 0. In analogy with (11) in that case we believe in the validity of the 
following conjecture. 
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Conjecture 2 For p — 0 the representing measure is given by 


W(x) = / e~ tx <! M ] 

i[log(l/p),oo) 


e ™ < /Wi og (i/p) + f e 4 


<?,p 


(B, b) 

(A, a) 


dt, 


where A is a positive constant and hi og (i/p) denotes the unit mass concentrated at the point 

log(l/p)- 

Unlike the case of (S'-function (11), no expansion similar to Nor lurid’s is available for H- 
function. We plan to study its behavior in the neighborhood of the singular point log(l/p) 
in our future work. 

Corollary 2 Suppose p > 0, a, b > 0 and conditions (12) and (16) are satisfied (in 
particular, (12) is true under Theorem 5). Then the function 


t —y 


nj=i r (h 3 

nun* 


■ H qS I e-* 


(B, b) 

(A, a) 


(24) 


represents an infinitely divisible probability density on (0,oo). 


Proof. The result follows from Theorem 7 combined with [20, Definition 5.6] and [20, 
Theorem 5.9]. See also [12, Theorem 7]. □ 

Remark 5. The case A = B is also connected to probability as follows. Suppose 
k — 1,.., ,p are independent beta-distributed random variables, so that ( k has the density 


r(«fc + fik) 


Then the random variable u = ]Qfc=i C k k with A k > 0, has the following moments 


E( u x ~ l ) — TT @ k 2 TT r (A k x + a k -A k ) 

J- J. TT-L J. 

k =1 


^ r(ofc) ^ T(A k x + a<k + fik — A k ) ’ 
The probability density of u is expressed via Fox’s //-function: 


x > 0. 


pj r(a fc + 0 t ) 

Hl’fi [ 2 

k =1 


r (a k ) 


p,p 


(Aj, otj + fij - Aj),j = 1 
(Aj, a.j — Aj),j = 1 


See [15, Section 4.2.1], 


Theorem 8 Suppose A.B>0, a, b>0, p>0,p<l and Y7i =i A = S"/=i Bj- Then 
for all x G (0, p) the following identity holds: 


Kp * * 


(B,b) 
(A, a) 


log (p/x) J x/p 


H P A - 


<up 


x 


u 


(B.b) ^ Q( 
(A, a) 


du, 


(25) 


where Q(u ) is defined in (14). 
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Proof. Since — log(p) = 5b — 5a, where 5a (5b) stands for the entropy of A (B), we adopt 
the notation AS = 5b — 5a = — log(p). Further, denote 


m = 



t > AS 
t < AS ’ 


and define the (signed) measure 


du(u) = A S5 0 + P(u)du, 

where <5 0 is a unit mass at zero. It follows from the proof of Theorem 4 and the asymptotic 
formula (19) that representation (17) is true under conditions of the theorem. In terms of 
the measure du (17) takes the form: 


W\s) = W{log W) 1 = -W e~ su du{u). 

J [0,oo) 

On the other hand, under the hypotheses of Theorem 6 we have the representation 
W{s)= f e~ ts H{e~ t )I{t)dt, where H(x) = Hf (x 

Jo V ( A ’ a ) 


Hence, 


W'(s) = 


te ts H{e l )I{t)dt. 


Putting these representations together yields 


te- ts H{e~ t )I{t)dt= / e~ ts H{e~ t )I{t)dt 


' [0,oo) 


e~ su dv{u). 


The convolution theorem for the Laplace transform [22, Theorem 11.4] then leads to the 
formula 

tH{e~ l )I{t) = / H{e^ ( ' t ^ T ' l )I{t — r)dv{r) for t G (A5, oo). 

■A o,0 

Recalling the definition of du we can rewrite the above relation as 


rt—As 


tH{e~ t ) = ASH{e ) + 


H{e~ t+T )P{r)dT 


or 


rf~AS 


H{e~ t ) = 


t-AS 


H{e~ t+T )P{T)dT. 


Substituting x = e t and u = e T yields identity (25). 

Theorem 4 shows that under additional restrictions (16) we have the implication 


Q(t) > 0 on [0,1) 


K P < * 


(B, b) 
(A, a) 


> 0 on (0, p). 


□ 


The following much stronger assertion is supported by numerical evidence: 
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Conjecture 3 Suppose A, B > 0, a, b > 0, p > 0, p < 1 and Yli=i A t = J2j=i Bj- 
Then 


# 


zeros of 



(B,b) 
(A, a) 


on (0, p) }> < # | zeros of Q(t) on (0,1) 


Using Theorem 8 we were able to demonstrate that this conjecture is true for the case when 
the left hand side is equal to one. 
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